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Spatial search by quantum walk 
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Grover's quantum search algorithm provides a way to speed up combinatorial search, but is not 
directly applicable to searching a physical database. Nevertheless, Aaronson and Ambainis showed 
that a database of N items laid out in d spatial dimensions can be searched in time of order \/A for 
d > 2, and in time of order %/~N poly (log A) for d = 2. We consider an alternative search algorithm 
based on a continuous time quantum walk on a graph. The case of the complete graph gives the 
continuous time search algorithm of Farhi and Gutmann, and other previously known results can 
be used to show that s/~N speedup can also be achieved on the hypercube. We show that full \f~N 
speedup can be achieved on a d-dimensional periodic lattice for d > 4. In d — 4, the quantum walk 
search algorithm takes time of order %/N poly (log A), and in d < 4, the algorithm does not provide 
substantial speedup. 



I. INTRODUCTION 

Grover's quantum search algorithm [1] is one of the 
main applications of quantum computation. Given a 
black box function f(x) : {1, . . . , N} — > {0, 1} satisfy- 
ing 
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x 7^ w 
x = w . 



(1) 



Grover's algorithm can find the value of w using of order 
\/A" queries, which is optimal [2]. On the other hand, no 
classical algorithm can do better than exhaustive search, 
which takes of order N queries. Therefore Grover's algo- 
rithm can be used to speed up brute force combinatorial 
search. It can also be used as a subroutine in a variety 
of other quantum algorithms. 

Grover's algorithm is sometimes described as a way to 
search an unsorted database of N items in time 0(yN). 
But the algorithm as originally proposed is not designed 
to search a physical database. Suppose we had N items 
stored in a d-dimensional physical space, and that these 
items could be explored in superposition by a quantum 
computer making local moves (a "quantum robot" [3]). 
Naively, it would seem that each step of the Grover al- 
gorithm should take time of order N 1 ^, since this is the 
time required to cross the database. Performing \ N iter- 
ations, we find that the search takes time of order N? + * , 
so no speedup is achieved in d — 2, and full speedup is 
achieved only in the limit of large d. 

However, it is possible to do better than this naive ap- 
proach suggests. In [4], Aaronson and Ambainis present a 
model of query complexity on graphs. Within this model, 
they give a recursive algorithm for the search problem 
that achieves full yJV speedup for a d > 3 dimensional 
lattice, and runs in time y/N log 2 N in d = 2. (It is 
obvious that no algorithm can get speedup in d = 1.) 
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In this paper we approach the spatial search problem 
using quantum walks. Since random walks are commonly 
used in classical algorithms, it is natural to consider a 
quantum analogue of a classical random walk as an algo- 
rithmic tool. Here we consider the continuous time quan- 
tum walk [5]. On certain graphs, this quantum walk can 
yield exponentially faster hitting times than its classical 
counterpart [5, 6]. Indeed, a recent result shows that the 
continuous time quantum walk can solve a certain black 
box problem exponentially faster than any classical algo- 
rithm [7]. 

Quantum walks provide a natural framework for the 
spatial search problem because the graph can be used 
to model the locality of the database. We present a 
simple quantum walk search algorithm that can be ap- 
plied to any graph. Our algorithm could be implemented 
within the model of [4], but is actually much simpler be- 
cause it uses no auxiliary storage space. For the case 
of the complete graph, the resulting algorithm is simply 
the continuous time search algorithm of Farhi and Gut- 
mann [8]. On the hypercube, previous results can be 
used to show that the algorithm also provides quadratic 
speedup [9, 10]. However, in both of these cases, the 
graph is highly connected. Here, we consider the case of 
a cZ-dimensional cubic periodic lattice, where d is fixed 
independent of N. We find full \/N speedup in d > 4 
and running time 0(y/N log 3 / 2 N) in d = 4. In d < 4, we 
find that quadratic speedup is impossible, so the contin- 
uous time quantum walk algorithm is never faster than 
the Aaronson- Ambainis algorithm. 

We note that it is also possible to construct a quan- 
tum analogue of a discrete time random walk [11, 12] 
(although the walk cannot take place directly on the ver- 
tices of the graph [13]). This type of walk has been used 
to construct a fast search algorithm on the hypercube 
[14], and more recently, on a d-dimensional lattice with 
d > 2 [15]. The latter result outperforms our continuous- 
time walk algorithm for d = 2,3,4. However, similar 
performance can be achieved by a modification of the 
continuous-time algorithm [16]. 

This paper is organized as follows. In Section II we 
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review the continuous time quantum walk and show how 
it can be used to approach the search problem. In Section 
III we review the results in the high dimensional cases 
(the complete graph and the hypercube), casting them 
in the language of continuous time quantum walks. In 
Section IV we present the results for finite dimensional 
lattices, and in Section V, we conclude with a discussion 
of our results. 



II. QUANTUM WALK 

The continuous time quantum walk on a graph is de- 
fined in direct analogy to a continuous time classical ran- 
dom walk [5] . Given an undirected graph G with iV ver- 
tices and no self loops, we define the adjacency matrix 

^ k ~ \ otherwise [Z> 

which describes the connectivity of G. In terms of this 
matrix, we can also define the Laplacian L = A — D, 
where D is the diagonal matrix with Djj = dcg(j), the 
degree of vertex j. The continuous time random walk on 
G is a Markov process with a fixed probability per unit 
time 7 of jumping to an adjacent vertex. In other words, 
the probability of jumping to any connected vertex in 
a time e is 76 (in the limit e — > 0). This walk can be 
described by the first-order, linear differential equation 

—^-=ll^L jk p k {t), (3) 

k 

where pj (t) is the probability of being at vertex j at time 
t. Since the columns of L sum to zero, probability is 
conserved. 

The continuous time quantum walk on a graph takes 
place in an TV-dimensional Hilbert space spanned by 
states \j), where j is a vertex in G. In terms of these 
basis states, we can write a general state \4>(t)) in terms 
of the N complex amplitudes qj(t) = (j\ip(t)). If the 
Hamiltonian is H, then the dynamics of the system are 
determined by the Schrodinger equation, 1 

k 

Note the similarity between (3) and (4). The continuous 
time quantum walk is defined by simply letting H = 
— 7L. 2 Then the only difference between (3) and (4) is 



1 We have chosen units in which h = 1. 

2 Here the sign is chosen so that the Hamiltonian is positive 
semidefinite. We have defined L = A — D so that for a lat- 
tice, L is a discrete approximation to the continuum operator 
V 2 . A free particle in the continuum has the positive semidefi- 
nite Hamiltonian H = —V 2 (in appropriate units). 



a factor of i, which nevertheless can result in radically 
different behavior. 

As an aside, we note that the Laplacian does not 
provide the only possible Hamiltonian for a quantum 
walk. Whereas (3) requires ^ Ljk = to be a valid 
probability-conserving classical Markov process, (4) re- 
quires H = to be a valid unitary quantum process. 
Therefore we could also choose, for example, H — — 7 A 
All of the graphs we consider in this paper are regular 
(i.e., deg(j) is independent of j), so these two choices 
give rise to the same quantum dynamics. However, for 
non-regular graphs the two choices will give different re- 
sults. 

To approach the Grover problem with a quantum walk, 
we need to modify the Hamiltonian so that the vertex w 
is special. Following [8], we introduce the oracle Hamil- 
tonian 3 

H w = -\w)(w\ (5) 

which has energy zero for all states except \w), which is 
the ground state, with energy — 1. Solving the Grover 
problem is equivalent to finding the ground state of this 
Hamiltonian. In this paper we assume that this Hamil- 
tonian is given, and we want to use it for as little time 
as possible to find the value of w. Note that this Hamil- 
tonian could be simulated in the circuit model using the 
standard Grover oracle 

U m \j) = (-l) 8 *»\j). (6) 

However, in this paper we focus on the continuous time 
description. 

To construct an algorithm with the locality of a partic- 
ular graph G, we consider the time-independent Hamil- 
tonian 

H = - 1 L + H w = -jL -\w)(w\ (7) 

where L is the Laplacian of G. We begin in a uniform 
superposition over all vertices of the graph, 

l-> = ^Eli>. W 

j 

and run the quantum walk for time T. We then mea- 
sure in the vertex basis. Our objective is to choose the 
parameter 7 so that the success probability \{w\ip(T))\ 2 
is as close to 1 as possible for as small a T as possible. 
Note that the coefficient of H w is held fixed at 1 to make 
the problem fair (e.g., so that evolution for time T could 
be simulated with 0(T) queries of the standard Grover 
oracle (6)). 



3 More precisely, we should use H w = — u;|w)(w| where ui is a fixed 
parameter with units of inverse time. However, we choose units 
in which ui = 1. In these units, 7 is a dimcnsionlcss parameter. 
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One might ask why we should expect this algorithm to 
give a substantial success probability for some values of 
7, T. We motivate this possibility in terms of the spec- 
trum of H. Note that regardless of the graph, \s) is the 
ground state of the Laplacian, with L\s) = 0. As 7 — > 00, 
the contribution of H w to H is negligible, so the ground 
state of H is close to \s). On the other hand, as 7 — > 0, 
the contribution of L to H disappears, so the ground 
state of H is close to \w). Furthermore, since \s) is nearly 
orthogonal to \w), degenerate perturbation theory shows 
that the first excited state of H will be close to \s) as 
7 — > for large N. We might expect that over some in- 
termediate range of 7, the ground state will switch from 
\w) to |s), and could have substantial overlap on both for 
a certain range of 7. If the first excited state also has 
substantial overlap on both \w) and \s) at such values of 
7, then the Hamiltonian will drive transitions between 
the two states, and thus will rotate the state from \s) 
to a state with substantial overlap with \w) in a time of 
order 1/(E\ — Eq), where Eq is the ground state energy 
and Ei is the first excited state energy. 

Indeed, we will see that this is a good description of 
the algorithm if the dimension of the graph is sufficiently 
high. The simplest example is the complete graph (the 
"analog analogue" of the Grover algorithm [8]) which can 
be thought of roughly as having dimension proportional 
to N. A similar picture holds for the (log 7V)-dimensional 
hypercube. When we consider a d-dimensional lattice 
with d independent of N, we will see that the state \s) 
still switches from ground state to first excited state at 
some critical value of 7. However, the \w) state does not 
have substantial overlap on the ground and first excited 
states unless d > 4, so the algorithm will not work for 
d < 4 (and d = 4 will be a marginal case) . 



III. HIGH DIMENSIONS 

In this section, we describe the quantum walk algo- 
rithm on "high dimensional" graphs, namely the com- 
plete graph and the hypercube. These cases have been 
analyzed in previous works [8-10]. Here, we reinterpret 
them as quantum walk algorithms, which provides moti- 
vation for the case of a lattice in d spatial dimensions. 



A. Complete graph 

Letting L be the Laplacian of the complete graph, we 
find exactly the continuous time search algorithm pro- 
posed in [8] . Adding a multiple of the identity matrix to 
the Laplacian gives 



L + NI = N\s)(s 
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FIG. 1: Energy gap and overlaps for the complete graph with 
N = 1024. 

Therefore we consider the Hamiltonian 

H = -jN\s)(s\-\w)(w\. (10) 

Since this Hamiltonian acts nontrivially only on a two- 
dimensional subspace, it is straightforward to compute 
its spectrum exactly for any value of 7. For 77V <C 1, 
the ground state is close to \w), and for jN >• 1, the 
ground state is close to \s). In fact, for large N, there is 
a sharp change in the ground state from \w) to \s) as 7 TV 
is varied from slightly less than 1 to slightly greater than 
1 . Correspondingly, the gap between the ground and first 
excited state energies is smallest for 7 A <~ 1, as shown 
in Figure 1. At -fN = 1, for N large, the eigenstates are 
-^(\w) ± \s)) (up to terms of order TV -1 / 2 ), with a gap 

of 2/y/N. Thus the walk rotates the state from |s) to \w) 
in time tz^/N /2. 

B. Hypercube 

Now consider the n-dimcnsional hypercube with N = 
2™ vertices. The vertices of the graph are labeled by n-bit 
strings, and two vertices are connected if and only if they 
differ in a single bit. Therefore the adjacency matrix can 
be written as 

(7) 

where cri is the Pauli sigma x operator on the jth bit. 

In this case, we again find a sharp transition in the 
eigenstates at a certain critical value of 7, as shown in 
Figure 2. The Hamiltonian can be analyzed using essen- 
tially the same method we will apply in the next section, 
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FIG. 2: Energy gap and overlaps for the hypercube with N = 
2 1() = 1024. 



together with facts about spin operators. The energy gap 
is analyzed in Section 4.2 of [9], and the energy eigen- 
states are analyzed in Appendix B of [10]. The critical 
value of 7 is 



1 n 

' On / j 



. - = - + 0(n 2 , 
2 n ^ \ r) r n 

r=l 



at which the energy gap is 



7N [1 + oin ~ 1)] 



(12) 



(13) 



and the ground and first excited states are ^=(|tu) ± |s)) 
up to terms of order 1 jn. Again, we find that after a time 
of order \/N, the probability of finding w is of order 1. 

IV. FINITE DIMENSIONS 

Having seen that the algorithm works in two cases 
where the dimension of the graph grows with N, we now 
consider the case of a d dimensional cubic periodic lat- 
tice, where d is fixed independent of N. The minimum 
gap and overlaps of \s), \ w) with the ground and first ex- 
cited states are shown in Figure 3 for d = 2,3,4,5 and 
N 1000. In all of these plots, there is a critical value of 
7 where the energy gap is a minimum, and in the vicinity 
of this value, the state \s) changes from being the first 
excited state to being the ground state. In large enough 
d, the \w) state changes from being the ground state to 
having large overlap on the first excited state in the same 
region of 7. However, for smaller d, the range of 7 over 
which the change occurs is wider, and the overlap of the 
I to) state on the lowest two eigenstates is smaller. Note 
that in all cases, \s) is supported almost entirely on the 



subspace of the two lowest energy states. Therefore, if 
the algorithm starting in the state |s) is to work at all, 
it must work essentially in a two dimensional subspace. 

In the rest of this section, we will make this picture 
quantitative. We begin with some general techniques 
for analyzing the spectrum of H using knowledge of the 
spectrum of the graph. We then show the existence of a 
phase transition in 7, and we show that for any d, the 
algorithm fails if 7 is not close to a certain critical value. 
Next we consider what happens when 7 is close to its 
critical value. In d > 4, we show that the algorithm gives 
a success probability of order 1 in time of order V^V, and 
in d = 4, we find a success probability of order 1/logiV 
in time of order y/NlogN. Finally, we investigate the 
critical point in d < 4 and show that the algorithm does 
not provide substantial speedup. 



A. Preliminaries 

In this section, we show how the spectrum of H can be 
understood in terms of the spectrum of L. An eigenvector 
of H, denoted \ip a ), with eigenvalue E a , satisfies 



H\i, a ) = {- 1 L-\w){w\)\i, a )=E a \il> a ) 



i.e., 



{- 1 L-E a )\i, a ) = \w){w\^ a ) 
The state \tp a ) is normalized, so |( , i/' |'i/' )| 2 

Ra = \{w\4> a )\ 2 

and choose the phase of \ip a ) so that 

(w#a) = V^- 
We wish to calculate the amplitude for success 

(w\e- im \s) =Y,Wi>a)(ipa\s)e- iE « t , 



(14) 

(15) 
1. Define 
(16) 



(17) 



(18) 



so we only need those \ip a ) with R a > 0. 

L is the Laplacian of a lattice in d dimensions, periodic 
in each direction with period TV 1 / 6 *, with a total of N ver- 
tices. Each vertex of the lattice corresponds to a basis 
state where x is a <i-component vector with compo- 
nents xj e {0,1,..., iV 1 ^ — 1}. The eigenvectors of — L 
are \<f>(k)) with 



<a#(fc)> 



ik-x 



where 

27rmj 

0,±1, 
0,±1, 



(19) 



(20) 



.,±i(JVVd-i) 
.,±!(ATV<*_2), 



N 1 ^ odd 



-7V 1 / d N 1 ^ 



even, 



(21) 
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nergy gap and overlaps for d-dimensional lattices 
1000. From top to bottom, d = 5, N = 4 5 = 1024; 

= 6 4 = 1296; d = 3, N = 10 3 = 1000; d = 2, 
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FIG. 4: The function F(E) for a d = 2 dimensional periodic 
lattice with iV = 16 vertices, at 7 = 1. 



and the corresponding eigenvalues are 
£(k) =2 (rf-^cos(fcj) 



(22) 



Since (</>(fc)|u>) ^ 0, from (15) we have 

(7£(*0 - E a){<t>(k)\lpa) + (23) 

for any fc. We can therefore rewrite (15), using (17), as 



lVv> 



(24) 



-lL - E a 1 

Consistency with (17) then gives the eigenvalue condition 

1 



H- 



w) = 1 . 



-7^ - E a 

Using (19), this can be expressed as 



(25) 



F{E a ) = 1 , ^) = i^^. (26) 

A typical function F(_E) is shown in Figure 4. This 
function has poles where E = r )£{k). For E 7^ j£(k), 
(26) shows that F'(E) > 0, so there is an eigenvalue of 
H between every adjacent pair of eigenvalues of — 7L. 
Since F(E) — > as E — > ±00, there is also one nega- 
tive eigenvalue of H (corresponding to the ground state). 
Note that in the case shown in Figure 4, the eigenvalues 
£ = 2, 4, 6 of — 7L have degeneracies 4, 6, 4 because of the 
symmetry of the lattice. It follows that there are 3,5,3 
eigenvectors of H with eigenvalues E a = 2,4,6, all with 
(u>|V>a) = and thus not relevant to our purpose. These 
11 eigenvectors, together with the 5 relevant ones, make 
up the necessary total of 16. 



G 



The normalization condition on \ip a ) gives 
1 



Ra(w\ 



(- 7 L - E a f 



\w) = l 



i.e. 



F'(Ea) 



(27) 



(28) 



We also need the overlap of \ip a ) with \s). Since L|s) = 0, 
from (24) we have 



so that 



IWa>| 2 



V fig 



(s|w) , 



NE%F'(E a ) ■ 



Using (18), (24), and (25), 



-iE a t 



At t = 0, this gives the sum rule 

y 1 = i. 



(29) 



(30) 



(31) 



(32) 



We will see that the spectrum of H depends signifi- 
cantly on the behavior of the sums 



«? - 1 V 1 



(33) 



If d > 2j, then Sj t d can be approximated by an integral 

4 



as 



where 



S it d = Ij,d + o(l) 



d d k 



(2nf J_ v [£{k)V ■ 



(34) 



(35) 



The condition d > 2j is necessary for I^d to converge at 
k = 0. The numerical values of I\ t d and 7 2 ,d for d < 10 
are given in Table I. Note that I^d can also be calculated 
using the formula [17] 



1 



da • 



(2dy J (7-1)! 



[^(a/d)]" 1 (36) 



d 


h,d 




3 


0.253 




4 


0.155 




5 


0.116 


0.0184 


6 


0.0931 


0.0105 


7 


0.0781 


0.00697 


8 


0.0674 


0.00504 


9 


0.0593 


0.00383 


10 


0.0530 


0.00301 



TABLE I: Numerical values of the convergent integrals. The 
result for ii,3 is given exactly in [18]; the rest were computed 
numerically. 



where X is a modified Bessel function of the first kind. 

On the other hand, if d < 2j, then Sj t d can be well 
approximated by the contribution from values of k small 
enough that £{k) is approximately 



(37) 



(where we have used the notation k 2 = k\ + ■ ■ ■ + k%). 
Then 



where 



= 1 y 1 



(38) 



(39) 



Here the sum is over all values of the (i-component vector 
of integers to other than to = 0, and converges for large 
to 2 . Numerically, we find 



c 2 2 = 0.00664, c 2 3 = 0.0265. 



(40) 



In the borderline case d = 2j, Ij t d diverges logarithmi- 
cally at k 2 small and Cj.d diverges logarithmically at to 2 
large. In this case 



S 



1 



j,2j 



(4tt)J jl 



-hiN + 0(l). 



We will need 



S 12 = —\nN + A + 0(N- 1 ) 
An 



'2.4 



32tt 2 



\nN + 0(l) 



(41) 

(42) 
(43) 



where A = 0.0488 (the case j = 1, d = 2 is treated in 
greater detail in [19]). 



4 The little-o notation /(TV) = o(g(N)) means 
linijv^oo f(N)/g(N) = 0. In contrast, the more familiar 
big-O notation f(N) = 0(g(N)) means there exist constants 
c,N such that for all N > N , \ f(N)\ < c\g(N)\. 



B. Phase transition 

In this section, we show that the overlap of the state \s) 
on the ground or first excited state of H exhibits a phase 
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transition at a critical value of 7 for any dimension d. In 
fact, away from the critical value, \s) is approximately an 
eigenstate of H, so Schrodinger evolution according to H 
does not change the state very much. In the next section, 
we will show that the algorithm indeed fails away from 
the critical value of 7, and in the following sections we 
will consider what happens near the critical point. 

For 7 larger than the critical value (which will be de- 
termined below), the ground state energy is very close to 
0. This can be seen as follows. The eigenvalue condition 
(26) for the ground state energy E , which is negative, 
gives 



1 = F(E ) 



1 



< 



N\E \ N^ ol £(k) + \E \ 
1 1^1 

1 h± 
N\E Q \ + 7 



(44) 
(45) 
(46) 



where in the last line we have assumed d > 2. In this 
case, for 7 > (which will turn out to be the critical 
value), up to small terms, 



I So I < 



1 



7 



7 - h,d ' 
Using (30), we have 

IWo)| 2 = \l + E 2 Q J2h£(k) + \E \y 



k^O 



> 



- 2 ^ [£{k)f 



k^O 



> 



7 2 ^ [S{k)f 



(47) 

(48) 

(49) 
(50) 



Inserting the behavior of S2, d from (33), (38), and (41) 
and using the bound (47), we find 



l-|Wo)| 2 < 



1 



(7 - h,df 



0{N- V ) d>4 
• < 0( AT" 1 log AT) d = A 
0{N- 2 / 3 ) d = 3. 

(51) 

This shows that if 7 = I\ t d + e for any e > 0, then 1 — 
K'SlV'o)! 2 approaches zero as N — > 00. 

If d = 2, then I\p is logarithmically divergent, but us- 
ing (42) in (45) we can apply a similar argument when- 
ever 7 > In N + A, in which case we have 



\Eo\ < 



7 



AT 7- ^lniV- A 



(52) 



and 



l-|Wo}| 2 < 



( 7 -^lnAT-A)2 



xO(l). (53) 



This shows that if 7 > (± + e)lnN 1 then 1- \(s\ip Q )\ 2 < 
l/(e\nN) 2 , which approaches zero as N — > 00. 

Similarly, for d > 2 and for 7 < the first excited 
state \ipi), with energy E x > 0, is essentially \s). Here 
we find 



'-f*'-iffii + s2?Fii (54 » 



> 



1 1^1 
1 



AT£i 7 ' 
so that, up to small terms, 

1 7 



£1 < 



A^ h. d - 7 
Again applying (30), we find 



(55) 
(56) 

(57) 



1-I(*>| 2 < 



{h,d - i f 



0{N- V ) d>4 
x \ OiN- 1 log N) d = 4 
0{N- 2 ' Z ) d = 3. 

(58) 

We see that 7 = l\.d is the critical point. In d = 2 
we can apply similar reasoning to obtain that for 7 < 
^lniV + A,' 



Wi>| 2 < 77 



1 



xO(l). (59) 



In this case 7 = In A' + ^4 is the critical point. 



C. Failure of the algorithm away from the critical 
point 

In this section we will show that the algorithm fails 
away from the critical point, regardless of dimension. 
The results (51) and (58) are actually sufficient to show 
that away from the critical point in d > 4, the algorithm 
can be no better than classical search, but we will give a 
different argument for consistency of presentation. 

First we consider the regime where 7 is larger than the 
critical value. In the previous section, we saw that in 
this case, the ground state energy E is small. This is 
sufficient to imply that the success probability is small 
at all times. Combining (31) and (32), we see that the 
amplitude at an arbitrary time must satisfy 



(w\e- im \s) \ < 



< 



(- 



1 

7A7 \\E \F'(Eo) 
2 



- 1 



(60) 
(61) 



VN\Eo\F'(E ) ' 
Furthermore it is clear from the definition of F(E) that 

1 



F'(E ) > 



NE 2 ' 



(62) 
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so 



\(w\e- lHt \s}\ <2VN\E a \. 
Using (47), we find that for d > 2, 

2 7 



<itf|e-""|a)| < 



7 ~ h,d ' 



(63) 



(64) 



This shows that if 7 = + e for any e > 0, the success 
probability is never more than a constant factor larger 
than its initial value, no matter how long we run the 
algorithm. If d = 2, then 7i j2 is logarithmically divergent, 
but using (52) we find 



(w\e- im \s)\< 



7 



(65) 



This shows that the algorithm fails if 7 > (j^ + e) In N 
for any e > 0. 

Now we consider the case where 7 is smaller than the 
critical value. For d > 4 and E < 0, we have 



F(E) 



1 



(27r) rf 7 7 £(fc) + |£| 
1 f d d k \E\ 



(66) 



d d k 



(2w) d J j£{k) {2n) d J 1 £{k)[ 1 £{k) + \E 
h. d \E\ f d d k 



> 7 7 2( 2 ^)rf j [f(fc)]2 
7 7 2 

Using the fact that F(E ) = 1, this shows that 

7<7m - 7) 



(67) 
(68) 

(69) 



\En\ > 



t2.d 



From (16) and (28), it is clear that F'(E) > 1, so using 
(61) gives 



(w\e- iHt \s) \ < 



21 



2,d 



VN j(h,d - 7) 



(71) 



A similar argument can be used for d = 3,4. With 
d = 4, we have 



F(E) 



1 



d A k 



(2tt) 4 J j£(k) + \E\ 
1 r d 4 fc |£| 



(72) 



(27r) 4 7 7 £(fc) (27r) 4 7 7 £(fc)[ 7 £:(fc) + |£ 



d 4 fc 



> 



h. 



\E\ 



2tt 



fcdfc 



7 32 7 7 lip + \ E \ 
7 256 7 2 n \ \E\ 



(73) 
(74) 

(75) 



where the third line follows because cosfc < 1 — 2(fc/7r) 2 
for |fc| < 7r, which implies £{k) > -^k 2 . We have also 
used the fact that k 2 < dn 2 to place an upper limit on 
the integral. This shows that for any e > (with e < 1), 
there exists a c > such that 



F(E) > 



h,4 c\E\ 



7 7 



2-e 



so that 

\Eo\ >c' 7 (/ M -7) 1/(1 - £) 
for some d > 0, and therefore 

1 2 



IHe"^>)|< 



With d = 3, we have 



d 3 fc 



7 £(fc) + l£l 

d 3 fc \E\ 



(76) 



(77) 



(78) 



(79) 



d 3 /c 



> 



(2tt) 3 7 7 £(fc) (2tt)3 J 7 £(fc)[ 7 £(fc) + 

(80) 

/1.3 i*i r d fc (gi) 

(82) 



7 87 7 ^ 2 + \E\ 



7 32 7 3 / 2 



where in the third line we have again used £ (k) > ^?k 2 . 
In this case we find 



hp 1 1024 ,t ^2 
l*o| > — J-7(-fl,3 - 7) 



(70) 

which shows that 



(s\e- im \w)\ < 



2tt 4 



ViV 1024 7 (7i i3 - 7 ) 2 



(83) 



(84) 



Finally, with d = 2 we use a different argument. Here 
we have 



F'{E) 



1 



d 2 fc 



> 



(2^) 2 7 [j£ (k) + \E\} 2 
kdk 



1 r 

2^ Jo 



2nJ ( 7 fc 2 + |E|) 2 

7T 

4|£|(|£| +^ 2 7 ) 



(85) 

(86) 
(87) 



where the second line follows since cos k > 1 — \ k 2 , which 
implies £{k) < k 2 . In the second line we have also used 
the fact that the entire disk |fc| < n is included in the 
region of integration. Equation (87) shows that 



\E\F'{E) > 



4(\E\ +7T 2 7 ) ' 



(88) 



9 



so that 



|He |s>| < 



1 8(|g |+^ 2 7 ) 

VN 7T 



(89) 



which is (9(1/ \/N) for 7 = 0(1), and 0({logX)/\/N) for 
any 7 < ^ In X + A. 

The arguments for the case where 7 is smaller than 
the critical value can be made tighter by a more refined 
analysis. For example, by considering the behavior of 
F'(E), one can give a bound whose dependence on I\.d— 7 
is linear for all d > 2, not just for > 4. Futhermore, 
the careful reader will note that our bounds for d > 2 all 
become useless as 7 — > 0, but it is easy to see that the 
algorithm cannot be successful for small values of 7. 

Altogether, we see that the algorithm cannot work any 
better than classical search if 7 is not chosen close to its 
critical value. It remains to investigate what happens 
near the critical point. 



D. The critical point in d > 4 

In this section we investigate the region of the critical 
point in the cases where the algorithm provides speedup. 
First we consider the case d > 4. Separating out the 
k = term in (26), we have 



f(e) = — L + ly . 



1 



XE 1 Nf^ ol £{k)-E 



(90) 



If \E\ <C 7<£(fc) for all k ^ 0, then for large N, we can 
Taylor expand the second term to obtain 



which gives 



F'(E) 



XE 2 



h.d 

7 2 



(91) 



(92) 



The critical point corresponds to the condition 7 = Ii t d- 
At this point, setting (91) equal to 1 gives two eigenval- 
ues, 



Eq 



h.d 



^fh~d~x 



+ 



(93) 



which correspond to the ground and first excited state, 
with a gap of order A" 1 / 2 . Since £{k) w (27r) 2 A- 2 / d 



for 



= 1, we see that the assumption Eq, E\ -C 7 £(&) 



holds for all k 7^ 0. Furthermore, for the ground and first 
excited states at 7 = Ii^, (92) gives 



F'(E ) « F'(Ei) 



2J24 

I 2 
1 i,d 



(94) 



the first excited state is small, since as can be seen using 
(32) we have 



— y 1 



1 



1 



1 



1 



E Q F'(E ) E 1 F'(E 1 



(95) 



Using (93) and (94), we see that the O(vJV) contribu- 
tions from 1/E F'(E ) and 1/E\F'{E\) cancel, so the 
right hand side of (95) is o(l). Thus, using (31), we find 



(w\e 



-iHt 




sin 



h,dt 
VhJX, 



(96) 



The success probability is of order 1 at t = I24X / I\ t d- 
Straightforward analysis shows that a similar condition 
holds so long as 7 = I\.d ± 0(A -1 / 2 ), exactly the width 
of the region that cannot be excluded based on the argu- 
ments of Section IV C. 

In d = 4, I2, d does not converge, so the result is modi- 
fied slightly. In this case (91) holds with l^.d replaced by 
— In A, so the ground and first excited state energies 
are given by 



Eq ft! — - 



h, 



^N\nX 
and we find 



Ei 



XlnX 



, (97) 



F'(E ) w F'(Ei) 



In A 

167T 2 J? >4 ' 



(98) 



Therefore 



\(w\e 



iHt 1 



IlA 



\nX 



sin 



'1,4 ' 



■XlnX 



(99) 



which shows that running for a time of order y/N log 
gives a success probability of order 1/ log TV. Using 
O(logA) repetitions to boost the success probability 
close to 1, we find a total run time 0{yfX log 3 / 2 TV). 5 
One can show that simila r conditions hold so long as 
7 = /i,4±0( v / (logA)/A). 

For d < 4, the expansion (91) fails to find states whose 
energies satisfy E <C j£{k). Indeed, we will see in the 
next section that the algorithm provides no substantial 
speedup in these cases. 



Now we want to use (31) to compute the time evolution 
of the algorithm. The contribution from all states above 



5 In fact, we could improve the run time of the algorithm to 
0(yN log AT) using amplitude amplification [20]. 
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E. The critical point in d < 4 

To handle the case d < 4, we rearrange the eigenvalue 
condition to extract the O(l) contribution to F(E): 



F<yE ^ NE + 1 £{k) 



-y 



E 



E\ 
(100) 

In d = 3, we can replace the middle term by h^/j 
for large N. To explore the neighborhood of the critical 
point in d = 3, we introduce rescaled variables a, x via 



7 : 

E : 



h,3 + 

47T 2 fi,; 

TV 2 / 3 



JVV3 



■ a; . 



(101) 
(102) 



Since the sum in the third term of (100) only gets sig- 
nificant contributions from small energies, we use (37) to 
give the approximation 



JV2/3 



(103) 



and we can analyze the sum using the same techniques 
we applied to calculate Sj^ in the case d < 2j. Then we 
have, for large N, 



F(E) w 1 



G 3 (x) 



/i, 3 JVV3 



where 



G 3 (x) = 



Air 2 



E 



to 2 (to 2 



(104) 



(105) 



Here the sum is over all integer values of to, as in (39), 
and similarly converges for large m 2 . The eigenvalue con- 
dition in terms of x is G 3 (x) = a, which has one nega- 
tive solution xq. Since G 3 (x) is independent of N, x 
is independent of N, and the ground state energy Eq is 
proportional to iV~ 2 / 3 . 

As we saw in Section IV C, a very small ground state 
energy implies that the success probability is small at all 
times. Using (63), we find 



|He |s>| < 



87r 2 Ji, 3 |a;o| 
ATi/6 



(106) 



Therefore the success probability is small no matter how 
long we run the algorithm. This fact is sufficient to im- 
ply that the algorithm cannot produce full square root 
speedup. Taking the time derivative of (31), we see that 



^IHe-^| S >| < 



which implies that 



ch He 



-iirt i 



t > I (io|e' 



-iHt\ 



i)\VN 



^ <107 » 



(108) 



Thus the time required to find w using classical repetition 
of the evolution for time t is of order 



t 



> 



N 



> 



]V 2 / 3 



87r 2 7 1>3 ko| 



(109) 
(110) 



regardless of t. In other words, the algorithm cannot 
produce full speedup. 

Similar considerations hold in the case d = 2. In this 
case, the critical point is at 7 = \nN + A, so we choose 



— lniV + A + a 
4ir 



2n\nN 

N 



7 : 

E = ., . 

In this case, we find 

F(E) w 1 



G 2 (x) - a 
±AnN 

4-rr 



(111) 

(112) 
(113) 



where G2(x) is defined as in (105), but with m having two 
components instead of three. Again we find a solution 
x < that is independent of N, and applying (63) gives 



Me-^| S )|< 47r|a;o|lniV 



Vn 



(114) 



(Note that we could have reached a similar conclusion 
using (89).) Using (109), we find 



t 



\{w\e- 



iHt\ 



> 



N 



4tt\x \ logN 



(115) 



so the algorithm also fails near the critical point in d = 2. 



V. DISCUSSION 

In this paper we have presented a general approach 
to the Grover problem using a continuous time quantum 
walk on a graph. We showed that quadratic speedup can 
be achieved if the graph is a lattice of sufficiently high 
dimension (d > 4). Although we had originally hoped 
to find a fast algorithm in d = 2, we found that our 
approach does not offer substantial speedup in this case. 

Our algorithm begins in the state |s), which is dclocal- 
ized over the entire graph. One might demand instead 
that we start at a particular vertex of the graph. How- 
ever, it is clear that \s) can be prepared from a localized 
state using 0{N 1 / d ) local operations. In fact, we could 
also prepare \s) by running the quantum walk search al- 
gorithm backward from a known localized state for the 
same amount of time it would take to find \w) starting 
from \s). 

The quantum walk search algorithm is related to a 
search algorithm using quantum computation by adia- 
batic evolution. Adiabatic quantum computation is a 
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way of solving minimization problems by keeping the 
quantum computer near the ground state of a time- 
varying Hamiltonian [9] . In the adiabatic version of the 
search algorithm, the quantum computer is prepared in 
the state \s) (the ground state of H with 7 large), and 7 
is slowly lowered from a large value to 0. If 7 is changed 
sufficiently slowly, then the adiabatic theorem ensures 
that the quantum computer ends up near the final ground 
state \w), thus solving the problem. The time required to 
achieve a success probability of order 1 is inversely pro- 
portional to the square of the gap between the ground 
and first excited state energies. On the complete graph, 
the fact that the gap is only small (of order N^ 1 / 2 ) for 
a narrow range of 7 (of order TV" 1 / 2 ) means that 7 can 
be changed in such a way that time 0(y/~N) is sufficient 
to solve the problem [21, 22]. Since the gap has similar 
behavior for the hypercube and for <i-dimcnsional lat- 
tices with d > 4, quadratic speedup can also be achieved 
adiabatically in these cases. In d = 4 t he gap is o f order 
1 / y/N \ogN for a range of 7 of order \J (log N)/N, so the 

run time is again 0(sfN 'log 3 / 2 N). In d < 4, no speedup 
can be achieved adiabatically. 

Yet another way to solve the Grover problem uses a 
sequence of measurements of H. For any adiabatic al- 
gorithm, there is a related algorithm that uses only a 
sequence of measurements to remain in the ground state 
of a slowly changing Hamiltonian [10]. The case of a 
hypercube was presented in [10], and our present results 
show that this algorithm can also be used when the graph 
is a lattice with d > 4. However, to realize the measure- 
ment dynamically, the Hamiltonian H must be coupled 



to a pointer variable, which must be represented using 
auxiliary space. 

Although the quantum walk algorithm does not per- 
form as well as the Aaronson-Ambainis algorithm in 
d = 2,3,4, it does have certain advantages. The quan- 
tum walk algorithm uses simple, time-independent dy- 
namics rather than a recursive procedure. Furthermore, 
the quantum walk algorithm uses only a single basis state 
for each vertex of the graph, whereas the algorithm of [4] 
needs substantial auxiliary space. 

The actual complexity of the search problem in d = 2 
remains an open question. It would be interesting ei- 
ther to improve on the algorithms of [15, 16] or to 
prove a lower bound showing that full speedup cannot 
be achieved. 
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